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PREFACE 


The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing Is a 6-year program of research, development, evaluation, and 
application of aerospace remote sensing for agricultural resources, which began 
in fiscal year 1980. This program Is a cooperative effort of the National 
Aeronautics and Space Administration, the U.S. Agency for International 
Development, and the U.S. Departments of Agriculture, Commerce, and the 
Interior. 

The work which is the subject of this document was performed within the Earth 
Resources Research Division, Space and Life Sciences Directorate, at the Lyndon 
B. Johnson Space Center, National Aeronautics and Space Administration. Under 
Contract NAS 9-15800, personnel of Lockheed Engineering and Management Services 
Company, Inc., performed the tasks which contributed to the completion of this 
research. 
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1. INTRODUCTION 


WATBAil Is a soil water point profile model developed by Dr. C. H. M. Van Bavel 
at Texas A&M University under contract to the NASA Johnson Space Center 
(ref. 1). This model Is patterned after the soil -pi ant-atmospheric model pre- 
sented earlier by Van Bavel and Ahmed (ref. 2). The computer program for the 
model Is coded In IBM's latest Continuous System Modeling Pregram III (CSMPIII, 
refs. 2 and 3). The use of CSMPIII In soil water modeling has been reported by 
Dr. Van Bavel and his colleagues In a number of papers (ref. 4). 

The CSMPIII has been developed primarily to solve the nonlinear partial 
differential equations of dynamic systems. The program has many calculation 
and printing capabilities In addition to standard FORTRAN routines. The 
general use of CSMPIII In soil water dynamics has been presented In a book by 
HI 1 lei (ref. 5) In which the basic soil water model and the description of the 
capabilities and the use of CSMPIII are well described. 

WATBAil has been designed to be general enough to represent realistically a 
wide range of soil -crop-atmospheric processes and conditions. In addition to 
the use of CSMPIII, WATBAil has several unique features that are not found In 
other soil water models. (See ref. 6 for the comparative characteristics of a 
number of soil water models.) For example, evaporation and transpiration are 
each determined separately and directly from the Input data. Also, the water 
flow through the crop Is determined by a difference In water potential divided 
by the appropriate crop resistance. Another feature is the determination of a 
canopy temperature from an energy balance approach. 

The purpose of this report Is to describe the model and Its output characteris- 
tics as the inputs are varied over realistic ranges. These characteristics are 
determined from a simulation study which can also indicate the boundaries for 
realistic simulation and the sensitivity of the output to given changes in the 
Input. A detailed summary of the model Is provided In the next section. 
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2 . MODEL SUMMARY 


2.1 GENERAL 

Any model of natural phenomena generally represents an approximation of the 
actual physical, chemical, or biological processes Involved. A discussion of 
the general processes involved in soil water changes are discussed in detail In 
reference 6. These processes can be represented by the following eouetloo: 

* P - I ♦ t>0 ♦ L - E - T ♦ F ♦ D fl) 

where 

aSW = the change in soil water for the layer for a given time interval 
P = precipitation or irrigation 
I 2 interception of P by the crop cover 

RO a surface net lateral flow, convergence (♦) (run-on), divergence (-) 
(runoff) 

L * net subsurface lateral movement, generally considered negligible 
E s evaporation (+) or condensation (-) 

T = transpi ration 

F * net vertical flux in layer, gain (+), loss (-) 

D * net flux at lower boundry, drainage (-), capillary rise (♦) 

Most of the terms on the right side, except precipitation, a^e not specified 
directly, but are estimated from functions involving atmospheric, pla"t, a f d 
soil parameters. WATBALl considers the processes represented by the following 
equation: 

:$W « o - E - T ♦ p * 0 '2) 

The actual representation of these processes is outlined in the schematic proc- 
ess flow diagram in figure 1. The boxes indicate processes modeled by fac- 
tions and logic steps. Input data are underlined; water losses are indicated 
by double arrows. The general flow of algor ’thm operations and data use is 
shown by the arrows between the boxes. Considerable feedback is indicated. 
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Figure 1.- Flow chart for VIATBAL1. 
















Soil water evaporation Is modeled by the processes Indicated on the left of the 
diagram. Transpiration and root uptake are modeled by the processes on the 
right. The middle section represents Intermediate and general calculations. 

The bottom section represents the processes Involved In determining the new 
soil water distribution. The calculations and operations represented here are 
discussed in detail below. A listing for the code for the model Is presented 
in Appendix A. 

The number of soil layers to be represented by the model and the thicknesses of 
the layers are arbitrary and can be varied to best represent the depths of 
Interest. The model can also group the layers Into larger units which have 
similar soil hydrologic properties. These properties _ moisture release and 
hydraulic conductivity data — are provided through tables of values. 

2.2 SOIL EVAPORATION 

The soil evaporation Is determined primarily by the net radiation at the ground 
under the canopy and the amount of water in the surface layer. The calculation 
is made by the following two equations: 

EVS = EV$1*EXP[PP0T(1)]/46.97*TAK (3) 

EVS1 = (EPS/EPSH)*NRBS/LH*1000.0 (4) 


where 

EVS = rate of evaporation from soil surface (m/s) 

EPS = ratio of a, which is the change in saturation vapor pressure with 

temperature, to y, the psychrometer constant 

EPSH = (EPS + 1) 

LH * latent heat of evaporation (J/kg) 

NRBS = net radiation at soil level (W/m^) 

PPOT(l) * matric potential of surface layer (m) 

TAK * temperature of the air (K) 
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PPOT(l) Is Initially an Input datum, but it Is then determined from the new 
water distribution values by the model calculations. TAK Is determined from 
the Input data on maximum and minimum Celsius temperatures by using linear 
Interpolation between the appropriate temperatures. Dewpoint values (DPTC) are 
determined In the same manner from the Input maximum and minimum dewpoints. 

EPS, LH, and NRBS are also computed from the appropriate Input data. The 
functions used are listed in table 1. NRBS is determined from a radiation 
energy balance calculation which is derived from the Incoming solar radiation, 
the long wave energy balance between air and canopy, the canopy albedo, and 
transmittance. 

2.3 TRANSPIRATION 

Transpiration is determined from the input data by the difference between the 
water vapor potentials of the atmosphere and the leaf, and the canopy resist- 
ance to the vapor flow. The following equation shows how this is determined: 

TRC * (HI - HA)/1000.0*RCW (5) 

where 

TRC * canopy transpiration rate (m/s) 

HL = absolute humidity of leaf interior (kg/m 3 ) 

HA * absolute humidity of the atmosphere (kg/m 3 ) 

RCW = total canopy resistance to water vapor diffusion (s/m) 

The total canopy resistance, RCW, is made up of two terms: Rt, the leaf 
stomatal resistance (s/m), and RA, the canopy resistance (s/m). The functions 
determining RL and RA, as well as HL and HA, are presented in table 2. As can 
be seen from table 2, the leaf stomatal resistance, RL, is determined by the 
current crop water potential, WPOTCR, from a table relating RL to WPOTCR. RA, 
the resistance to canopy diffusion, is a function of the daily mean wlndspeed 
which is an input datum. 

To determine a new WPOTCR, the following relationship Is used. (See table 2 
for definitions.) 
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TABLE 1.- FUNCTIONS FOR CALCULATING EVAPORATION PARAMETERS 
IN EQUATIONS (3) AND (4) 


LH = 2.49463 x 10 6 - 2.247 x 10 3 TL; latent heat of evaporation (.i/kg) 

EPS = 0.921 - 0.00262TL ♦ 0.00308TL 2 , where TL, Interim canopy temperature. Is 
calculated by an Implicit CSMPIII routine; l.e., TL * IMPL (TAC, 0.01, 
FTL) where TAC Is the air temperature determined from a linear Inter 
polatlon between the Input data on maximum and minimum temperatures. 

FTL = TAC - SHCA*(RA/SH); final canopy temperature # C 

RA * canopy resistance (see table 2) 

SH = 350089. 17/TAK; specific heat of air at constant volune (J/m 3 ) 

TAK = TAC + 273.16° K; see TAC above 

SHCA = LTR - NRBC; sensible heat transfer between canopy and atmosphere 

NRBC = (GRMABSC) + (1.0 - FTSR)*(SKL - LWRC); 

net radiation absorbed by the canopy (Vl/nr) 

LTR = (HL - HA)*LH/RCW; latent heat of transpiration; see table 2 for HL, HA, 
and RCU 

GR = (436.33 DGR/DL )*si ne[ ( STIME - 12. + DL/2)*ir/DL] 

DGR = dally total global radiation (mJ/m 2 ); dally Input data 
DL = day length (hours); dally Input data 
STIME = time of day (hours) 

ABSC = 0.0032 + 0.3084(LAI) - 0.05323(LAI 2 ) + 0.003667(LAI) 3 ; canopy 
absorptance 

LAI = leaf area Index; daily plant Input data 

FTSR = 0.9842 - 0.6755(LAI) + 0.1595(LAI) 2 - 0.0124(LAI) 3 ; view factor of 
diffuse radiation through the canopy 

SKL = o(TAK) 4 *(0.605 + 0.039 /I41TT.W) 

= long wave radiation from sky (W/m 2 ) 

HA = atmospheric absolute humidity (kg/m 3 ); see table (2) 
o = Stephan - Boltzmann's Constant 

LWRC « o(TL + 273. 16) 4 ; long wave radiation from canopy (W/m 2 ) 

NRBS = GR* (1.0 - ALBC - ABSC) + FTSR* (SKL - LWRC); 

= net radiation at soil level (W/m 2 ) 

ALBC = 0.124 - 0.009988 (LAI) + 0.007142(LAI) 2 - 0.000583(LAI) 3 , canopy albedo 
♦indicates multiplication 
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TABLE 2.- FUNCTIONS FOR CALCULATING TRANSPIRATION 


ML * 1.323 exp[17.27TL/(237.3 ♦ TL)]/(Tl + 273.16); absolute humidity of 
leaf 

HA - 1.323 exp[17.27DPTC/ (237.3 + OPTC)]/(273.16 + DPTC); absolute 
humidity of atmosphere 

RCVf * RL + RA; total canopy resistance (s/m) 

RA * AL0G(2.0/zq) 2 /0.16*SA; canopy resistance to diffusion 

SA = average dally wlndspeed (m/s); Input 

RL a i.o/RL'*(LAI); leaf stomatal resistance 

RL' Is obtained from an Input table relating RL to VIPOTCR 

WPOTCR = WPSEFF + WPCRMN - TRC*SRCR/LAI ; crop water potential (m) 

WPCRMN = crop water potential at zero transpiration and zero soil water 

potential (m); input 

SRCR a specific resistance to water uptake (s) 
m 

WPSEFF * £ PPOT ( I ) *RF ( I ) ; effective soil water potential (m) 

PPOT(I) = soil water matrlc potential In layer I (m) 

RF(I) = 2.0[^j- DEPTH(I)/RD 2 ]*TC0M(I); where RF( I) < 0, RF(I) = 0, root 
distribution parameter 
DEPTH(I) = depth of layer I; Input 
TCOM(I) = thickness layer I; Input 
RD a dally root depth (m); input 



( WPOTCR - MPCRMN) - WPSEFF * TRC*^B. 


( 6 ) 


This equation relates an effective potential difference on the left side of the 
equation to the product of the current water flow term (TRC) and a normalized 
resistance to water uptake. This equation Is used to find a new WPOTCR from the 
current transpiration and the current effective average soil water potential 
weighted by the root distribution. This latter parameter Is determined by the 
following equation. 

WPSEFF = £ PP0T(I)*RF(I) (7) 

1=1 


2.4 ROOT WATER UPTAKE 

The root water uptake, RC, Is now calculated by using the following equation: 

RC ( I ) = [(WPOTCR - WPCRN) - PPOT(I)]*UAI/SRCR)*RF(I) (8) 

This equation relates the water uptake from the layer I to a difference in 
potential which is divided by the specific resistance and weighted by the 
fraction of total roots in the layer. The difference in potential involves the 
effective crop potential minus the soil water potential of the layer. The root 
water uptake is related to the transpiration through the WPOTCR term as given 
by equation (6). 

2.5 RAINFALL AND INFILTRATION 

Increases in soil water are determined by the amount of the rainfall that 
infiltrates the soil. This infiltration is determined partly by the intensity 
of the rain and partly by the water amount of the near surface and surface 
layers. 

The intensity of the rain at a given time is regulated by the input data: time 

of beginning, time of ending, and rainfall total for the period. The program 
determines the midpoint of the period, distributes the rain linearly from zero 
at the beginning up to a maximum at the midpoint, and then distributes the rain 
linearly to zero at the end of the period. 
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The Infiltration Is determined by the following equations: 

OETAIN » INTGRKO.O, Rain • INFILT) (9) 

INCAP * [0.0 - HPOT ( 1 ) ]*0 • 5* [SATCON + C0ND(1)]/DIST(1) (10) 

where 

INTGRL = CMSPIII Integration function 
INCAP * maximum rate for Infiltration (m/s) 

HP0T(1) a water potential of surface lqyer (m) 

COND(l) = hydraulic conductivity of surface layer (m/s) 

SATCON * saturated conductivity (m/s) 

DIST(l) = distance from surface to midpoint of surface layer (m) 

OETAIN = amount of rainfall at surface that has not Infiltrated the soil 
INFILT = amount of rainfall that has Infiltrated the soil 

The net result of the equations and subsequent logic is that at a given time 
the Infiltration Is limited to DETAIN but all the rainfall eventually Infil- 
trates the soil. If the rate is greater than INCAP, the Infiltration Is spread 
over a longer time period. 

2.6 SOIL WATER PROFILE CHANGE 

The total amount of water in the profile or a layer during each time Interval 
is increased by infiltration, and it is decreased by evaporation at the 
surface, root uptake below the surface, and drainage at the lower boundary. 

The final step in the time interval is the determination of the net flux In 
each layer by the following equations (- = upward movement; + = downward 


movement): 

HLUX(I) = [HPOT(I-l) - HP0T(I)]*AVC0ND(1)/DIST(I) (11) 

FLUX(NLL) * COND(NL) (12) 

FLUX(l) = INFILT - EVS (13) 

NFLUX(I) = FLUX(I) - FLUX(I ♦ 1) + RC(I) (14) 
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where 

FLUX(I) = flux across the top boundary of layer I (m/s) 

FLUX (I ♦ 1) * flux across the bottom of layer I (m/s) 

HP0T(1) ■ total soil water potential (head) of layer I (m) 

AVCOND(I) = average hydraulic conductivity at boundary I (m/s) 

DISK I) * distance between midpoint layer (I - 1) and layer I (m) 

COND(NL) * hydraulic conductivity of last layer (m/s) 

FLUX(l) * flux at upper boundary 

INFILT * infiltration (m/s) 

EVS - evaporation (m/s) 

NFLUX(I) = net flux In layer I (m/s) 

RC(I) = water uptake from layer I (m/s) 

These equations show that the surface flux In a given time Interval Is the 
Infiltration minus the evaporation. The flux at the lower boundary (drainage), 
as determined by the program algorithm. Is equal to the conductivity for the 
layer and is always downward. The net flux In a layer Is equal to the differ- 
ences in the boundary fluxes minus the water uptake [-RC(I) = water loss]. To 
get the new soil water profile, the net flux in each layer I Is multiplied by 
the time Interval (seconds) and added to the contents of layer I at the 
beginning of the Interval. 

2.7 WATER BALANCE COMPUTATIONS 

As a check on the many calculations and operations in the model, a net balance 
value between the initial water amount in the profile, the resulting Infiltra- 
tion, the evapotransplration, the drainage, and the final water amount Is 
obtained using the following relation: 

BALANS * CUMWTR - ( IWATER + CUMINF - CUMETR - CUMDRN) (15) 
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where 

1WATER * the Initial amount of water In the profile (m) 

CUMWTR 3 new water amount In the profile after a given period of time (m) 

CUMINF 3 total amount of infiltration for a given time period (m) 

CUMETR 3 total amount of evapotranspl ration for a given time Interval (m) 

CUMDRN 3 total amount of drainage In a given time period (m) 

In order to appreciate the significance of the BALANS term. It Is necessary to 
analyze the term CUMWTR. This latter term Is the water amount In the profile 
at the beginning of the calculation plus the summation of the fluxes In the 
layers over the time period. This net flux Includes root uptake [-RC(I)j from 
each layer plus infiltration and evaporation at the upper boundary and d-ainage 
at the lower boundary. Relating this definition of CUMWTR to the terms ' n the 
BALANS equation above indicates that the BALANS term essentially compares the 
root uptake in the profile over a time interval (CUMRC) to the transpiration 
over the same time interval (CUMTR). Also included in the BALANS term a^e 
computational uncertainties resulting from the Initial computations and the 
ensuing integrations. If BALANS Is a positive number, CUMTR is generally 
larger than CUMRC; if negative, CUMTR should be smaller than CUMRC. In order 
to analyze the BALANS values, SUMRC and CUMRC were later added to the program 
code, and CUMRC was printed along with the other output. These parameters are 
similar to the other SUM and CUM parameters. 
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3. VALIDATION AND SENSITIVITY ANALYSES 


WATBAL1 Is a new, comprehensive model. The model should be tested extensively 
with field data so it can be used objectively, and the results can be Inter* 
preted with a known degree of confidence. However, if the needed field data 
are not available, which is the case here, a preliminary evaluation can be per* 
formed using simulated data. The use of simulated data can provide information 
on how well the model represents the generally anticipated characteristics of 
the domain modeled. 

If the simulated data are changed In a systematic manner, the variation In the 
output when compared to the variation In the Inout will also provide Insight 
Into the sensitivity of the output to uncertainties in the Input data. These 
sensitivity analyses can also indicate the accuracy and precision needed In the 
input data to obtain the desired accuracy and precision in the output. The 
data simulated represents the atmospheric and soil properties discussed below. 

3.1 ENVIRONMENTAL MODEL 

In order to perform the validation and sensitivity analyses, a standard data 
test set Is needed that represents typical conditions. In this regard the 
atmospheric and plant data for the standard data set are similar to the values 
used by Or. Van Bavel. The exception is that precipitation was not included in 
the standard set. The values are listed in tables 3 and 4. 

The soil characteristics are, however, different from those considered by Dr. 
Van Bavel. Basically, the properties of the Keith silt-loam profile near 
Colby, Kansas, were modeled. The layers were separated into two groups with 
different hydrologic properties in each group. The data for the hydrologic 
variable, moisture retention, are derived from three models: (1) regression 

(ref. 7), (2) Rogowski's (ref. 8), and (3) Ghosh's model (ref. 9). The basic 
input data for these models were obtained from previous soil surveys (ref. 10) 
or later in-situ measurements. The hydraulic conductivity data were then 
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TABLE 3.- STANDARD DATA SET 


Weather Input data 


JNM 

INPUT 

DL 

DGR 

TMAX 

TMIN 

OMAX 

DMIN 

SA 

BEGIN 

END 

RFT 

LAI 

RD 

121 . 

14.4 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.0 

0.5 

122. 

14.4 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.1 

0.5 

123. 

14.5 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.2 

0.6 

124. 

14.5 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.3 

0.6 

125. 

14.5 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.4 

0.7 

126. 

14.5 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.4 

0.7 

127. 

14.5 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.5 

0.7 

128. 

14.6 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.6 

0.8 

129. 

14.6 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.7 

0.8 

130. 

where 

JNM 

14.6 
* day 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

0.0 

0.0 

3.8 

0.9 


DL = day length (hours) 

DGR = daily global radiation (mj/m?) 


TMAX = maximum centigrade temperature 
TMIN * mlnlmjn centigrade temperature 
DMAX * maximum dewpoint centigrade temperature 
DMIN = mlnimun dewpoint centigrade temperature 
SA = mean dally windspeed (m/sec) 

BEGIN = beginning of rainfall (hour) 

END = end of rainfall (hour) 

RFT = amount of rainfall (m) 

LAI =leaf area Index 
RD * roofing depth (m) 
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TABLE 4.- PLANT AND PHYSICAL INPUT DATA 
(a) Parameters 


Sigma 

8 

5.67 X 10“ 8 

SATCON 

8 

0.30 X 10“ 6 

SRCR 

8 

1.00 X 10 9 

WPCRMN 

8 

-5.0 

WPOTCR (Initial) 

8 

-5.0 

FUNCTION RLVSWP 

8 

(0.0. 0.2), (50.0, .02), (150.0, 0.0002), 



(500.0, .0002), (20000.0, 0.000002) 


(b) Layer thickness and depth at midpoint 


Layer 

no. 

Thickness (m) 

Depth (m) 

1. 

1 1 

IjTF 1 1 

2. 



3. 



4. 



5. 


0.1143 

6. 

- ® | « v '. 

0.1524 

7. 



8 . 


gHHWfc&iH 

9. 

0.1524 

0.4572 

10. 

0. l c ?4 

0.6096 

11. 

0.1524 

0.7620 

12. 

0.1524 

0.9144 

13. 

0.1524 

1.0668 

14. 

0.1524 

1.2192 

15. 

0.1524 

1.3716 

16. 

0.1524 

1.5240 

17. 

0.1524 

1.6/&4 

18. 

0.1524 

1.8288 
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TABLE 5.- REGRESSION MODEL 
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• ••• 
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TABLE 7.- ROGOWSKI MODEL 
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obtained from Jackson's method using moisture retention data as Input 
(ref. 11). The moisture retention and hydraulic conductivity values derived 
for the three models are presented In tables 5, 6, and 7. 

Each simulation run was for a period of 10 days. Ten days were used mainly as 
a convenience since most of the runs In references 4 and 5 were for 10 days. 

In addition, a 10-day period Is often the Interval between soil moisture pro- 
file measurements. Simulation runs were made using the three models of the 
moisture retention and the standard atmospheric and plant data set described 
above. For each moisture retention model, simulation runs were then made using 
systematic variations from these standard conditions for a range of constant 
soil water values with depth (i.e., profiles). These constant profiles 
generally varied from 0.4 to 0.15 cm^/cm^ in 0.05 cm^/cm^ Increments. A large 
range of variation for most of the variables was used in conjunction with the 
regression model; fewer were used for the Rogowski and Ghosh models. Both crop 
and fallow conditions were simulated. 

The results of the simulations discussed above are presented below. The sensi- 
tivity characteristic for crops using the regression model will be presented 
and discussed first, followed by the Rogowski and Ghosh models. Then the 
responses of the model for fallow conditions are presented. 

3.2 CROP SIMULATIONS 

3.2.1 REGRESSION MODEL FOR WATER RETENTION 

3.2. 1.1 Daily E v , and T, and ET for Standard Conditions 

Figures 2, 3, and 4 present the simulated daily values of evaporation (E v ), 
Transpiration (T), and evapotranspiration (ET) over the 10-day period for wet, 
intermediate, and dry soils. These figures reproduce the observed three stages 
of drying, which are the constant rate stages (wet and dry) and the falling 
rate stage (intermediate). However, it can be seen that this three-stage 
drying is only present in the wet regime. In the other two regimes, only the 
falling rate and final constant rate stages are indicated. On the other hand, 
other simulat ons with lower solar radiation values (data not shown) extend the 
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initial constant rate stage to lower values of the initial soil water profiles. 
The cumulative evapotransplration for 10 days compared to initial soil water 
profiles are shown in figure 5. 
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and ET for 10 days for Intermediate soils. 
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Figure 4.- Daily values of E , T, and ET for 10 days for dry soils. 



0 days versus initial soil water profiles. 


The 10-day totals of ET, T, E v , and water loss ( IWATER - CUMWTR) for different 
Initial constant soil water profiles are presented In table 8 and figure 6. 

Total water loss Is higher than the ET values In the wet regime because of high 
drainage values. As the 10-day drainage decreases when e becomes smaller, the 
water loss curve becomes nearly Identical to the ET curve. It departs from the 
ET curve at still lower 0 values, and It Indicates less water loss than would 
be expected from the ET and drainage values. The amount of departure coincides 
with an Increase In the BALANS value since the model only allows drainage out 
the bottom boundary. 

An investigation of the cause of these high BALANS values has indicated that 
they are involved with the root uptake [RC(1)1 of soil water at Intermediate 
and low soil water values. The probable cause of this is discussed ir 9 later 
section. The response of E v , T, and ET to changes in the dally global 
radiation amount, D6R, Is shown In table 3 and figure 7. These results 
indicate that the model provides a nearly linear response to solar radiation 
changes over the range studied for the wet and dry boundary regimes. In the 
intermediate zone, the response is nonlinear. 

The curves indicate that for 0 = .40, a 10 percent change in the solar radia- 
tion at DGR = 20 provides approximately a 7.5 percent change in the 10-day ET. 
For 0 = .15, a 10 percent change in DGR gives approximately a 10 percent change 
in ET for 10 days. At a value of DGR = 10, the ET for 10-day response is 
approximately 9 percent for a 10 percent change in DGR at 0 * .40. For 0 * .15, 
the response is about 11 percent for a 10 percent change in DGR. 

The evapotranspiration on day 10 compared to the profile amount and daily solar 
radiation value is illustrated in figure 8. These curves show a consistent 
modified step character. 

The variation of ET for 10 days as the daily maximum temperature changes is 
presented in figure 9. The curves indicate that the response is nearly linear 
in the wet and dry regimes, but it is somewhat nonlinear in the intermediate 
regime. In general, the percentage response is less than the DGR response. 
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TABLE 8.- VAN BAVEL-CR OP-REGRESSION STANDARD 10-DAY TEST RESULTS 


Cumulative centimeters for 10 days 


Drainage BALANS 



Centimeters for last day 



Drainage 


0.40 

7.43 

5.20 

2.20 

0.901 

8.291 

0.036 

0.731 

0.536 

0.191 

.35 

7.14 

4.91 

2.21 

0.064 

7.141 

.059 

.666 

.469 

.194 

.30 

5.23 

2.93 

2.31 

0.004 

5.106 

.132 

.321 

.113 

.210 

.25 

3.39 

.96 

2.41 

3*10" 4 

3.228 

.166 

.285 

.073 

.212 

.20 

3.02 

.63 

2.39 

2xl0‘ 5 

2.610 

.410 

.278 

.073 

.205 

.15 

2.83 

.61 

2.21 

2xl0" 6 

1.586 

1.23 

.257 

.069 

.188 


0.40 

5.33 

3.94 

.35 

5.30 

3.91 

.30 

4.76 

3.36 

.25 

2.43 

.92 

.20 

1.93 

.44 

.15 

1.82 

.42 



0.906 

6.213 

0.026 

0.532 

0.410 


0.064 

5.325 

.039 

.522 

.401 


0.004 

4.681 

.080 

.394 

.271 

.125 

3xl0' 4 

2.301 

.129 

.183 

.050 

.132 

2xl0* 5 

1.606 

.328 

.179 

.050 

.130 

2xl0" 6 

.778 

1.04 

.170 

.048 

.120 





0.40 

3.35 

2.72 

0.64 

0.913 

4.239 

0.021 

0.337 

0.281 


.35 

3.36 

2.73 

.64 

.064 

3.402 

.019 

.337 

.281 

■si 

.30 

3.28 

2.66 

.64 

.004 

3.242 

.039 

.320 

.266 

.055 

.25 

1.80 

1.13 

.66 

3xl0“ 4 

1.678 

.121 

.106 

.047 

.059 

.20 

.95 

.28 

.67 

2*10‘ 5 

.616 

.272 

.089 

.031 

.059 

.15 

.90 

.27 

.64 

2xl0‘ 6 

.017 

.887 

.085 

.030 

.055 


II 


.64 
.64 
.64 
.25 1.37 


.20 .18 
.15 .18 



0.920 

3.255 

0.010 

HSI 

.064 

1.694 

.011 

.031 

.004 

3xl0" 4 

1.617 

.027 

.030 

1.280 

.086 

.025 

2xl0' 5 

-.05 

.233 

.024 

2xl0* 6 

-.611 

.789 


innaim* 







0.003 
.003 
.163 .003 

.119 .0025 

.017 .0022 

.017 .0021 
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Figure 7.- The response of ET for 10 days versus the dally global 
radiation amount (DGP.) and initial soil water profile. 
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Figure 9.- The variation of ET for 10 days versus the maximum 
temperature and initial soil water profile. 
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The response for changes In minimum temperature Is similar. The effects of 
changing the dewpoint are presented In figure 10. This simulation was accom- 
plished by lowering the minimum dewpoint. Changes In dewpoint result In con- 
siderably less response of ET for 10 days than changes in the temperature. 

The ET for 10-day response to dally mean wlndspeed changes Is shown In 
figure 11. These curves Indicate a nonlinear response for wet conditions and 
very little response for dry conditions. The decrease of ET for 10 days with 
Increase of wlndspeed when 0 » .15 Is unexpected and may reflect the RC Incon- 
sistency. Over most of the range of wlndspeeds, a 10 percent change In wind 
speed Indicates a 4 to 5 percent change In ET for 10 days for the wet boundary. 
Oryer soil conditions provide less of a change. 

The simulation response obtained for ET for 10 days from varying the LAI and RD 
In unison as a percentage of the standard values Is depicted In figure 12. 

These curves indicate that from small values of LAI and RD (l.e., shortly after 
f-uiergence to an LAI/RD of 15 percent) a 10 percent Increase In LAI/RD gives 
about c 13 percent change in ET for 10 days for the wet regime. The percent 
response in the drier regimes Is progressively larger. For progressively 
larger values of LAI/RD, the percent change In the ET response for a given 
change in LAI/RD progressively decreases in the wet regime. In the dry regime, 
the ET for 10 days actually decreases with further LAI/RD Increases, another 
unexpected response. 

The variation in ET for 10 days with the variation In the plant constant SRCR 
(the specific resistance to water uptake) Is Illustrated in figure 13. These 
curves indicate that the ET response to a given change in SRCR Is very small. 

In addition, other simulations Indicate that if SRCR Is decreased sufficiently 
at a given soil water amount the ET for the 10-day value changes sign, a very 
unrealistic situation (not shown on figure 3). Furthermore, Increasing the 
SRCR value above standard reduces the BALANS values for the drier regimes. 

Figures 14 and 15 show the response of two other plant constants, RLVSWP and 
WPCRM. The curves indicate negligible changes In the model reponse. The soil 
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Figure 10.- The response of ET for 10 days versus the dewpoint 
and initial soil water profile. 
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Figure 12.- The response of ET for 10 days from varying the LAI and RD in unison 
as a percentage of the standard values versus initial soil water profile. 



6 



Figure 13.- The variation in ET for 10 days versus constant specific 
resistance to water uptake (SRCR) and initial soil water profile. 
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Rl in RLVSWP Function 


Figure 14.- The response of ET for 10 days versus a change in the RLVSMP 
function as a percentage of standard values and intiial soil water 
profile. 
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water depth profiles on the 10th day are shown In figure 16. The wet profiles 
look realistic but become Increasingly unrealistic toward the dry regime. The 
profiles are unrealistic because the near-surface water values appear too high. 

3. 2. 1.2 Output Values Versus Hydraulic Conductivity 

Errors can occur In the hydraulic conductivity (h.c.) values since they are not 
measured, but calculated from the moisture release data by Jackson's method. 

To test responses to these errors, all the h.c. values were increased and 
decreased by 20 percent. The results of simulations using these changes, but 
otherwise standard Inputs, Indicate negligible effects on E v and T, but a 
definite effect Is noted on drainage and .*ater loss In the wet regime. 

3.2.2 R0G0WSK1 AND GHOSH MODELS FOR WATER RETENTION 

The results of simulations for a variety of conditions using Rogowskl's model 
are presented In table 9 and figures 17 and 18, and the results for Ghosh's 
model are shown In table 10 and figures 19 and 20. Inspection of these curves 
allows the following comments. 

The dally values for E v show very little difference from the regression model 
under similar environmental conditions. T values, however, are similar In the 
wet and dry regime, but remain high for a few more days before falling In the 
intermediate regime. 

The E y for 10 days as a function of 9 changes very little between the models 
for any of the environmental conditions simulated. On the other hand, the T 
changes are very small in the wet and dry regimes, but are considerable in the 
intermediate regimes; in some cases, changes of at least 200 percent occur. The 
shape of the ET curves, however, are similar to the regression model curves. 

Other significant differences in the output provided by the different water 
retention models occur in the values of the drainage and total water loss. As 
can be seen from table 8, 9, and 10, the regression model allows the least 
drainage and water loss, while the Rogowski model allows the most drainage and 
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TABLE 9.- VAN BAVEL-CROP-GHOSH STANDARD 10-DAY TEST RESULTS 



Cumulative centimeters for 10 days 

Centimeters for last day 

°v 

ET 

B 

m 

Drainage 

Water 

loss 

BALANS 

ET 

D 

B 

Drainage 


DGR = 30 


0.40 

7.38 

5.19 

2.19 

2.63 

9.907 

0.104 

0.734 

0.543 

0.191 


.35 

7.39 

5.19 

2.20 

.668 

8.044 

.017 

.742 

.547 

.192 

| 

.30 

7.23 

5.01 

2.21 

.084 

7.230 

.088 

.694 

.497 

.193 


.25 

4.55 

2.16 

2.36 

.006 

4.410 

.150 

.300 

.079 

.214 

iilSflSfl 

.20 

3.08 

.659 

2.42 

2*10" 4 

2.838 

.243 

.282 

.073 

.209 


.15 

2.84 

.611 

2.23 

4*10“ 6 

1.745 

1.10 

.257 

.069 

.188 

4xl0‘ 7 


DGR = 20 


0.40 

5.26 

3.88 

1.38 

2.68 

7.922 

0.017 

0.525 

1 0.405 


0.109 

.35 

5.28 

3.89 

1.38 

.670 

5.965 

-.020 

.529 

1 .409 



.30 

5.30 

3.91 

1.39 

.084 

5.323 

.056 

.526 




.25 

4.21 

2.74 

1.43 

.006 

4.088 

.130 

.321 

.184 


6xl0“ 4 

.20 

1.98 

0.46 

1.52 

2 * 10" 4 

1.794 

.182 

.182 


BPS 

2xlO" 5 

.15 

1.85 

0.43 

1.43 

4xl0" 6 

.991 

.857 


.048 

.122 



OGR = 10 


0.40 

3.27 

2.64 

0.634 

2.73 

6.014 


0.328 

0.273 



.35 

3.29 

2.66 

.635 

.672 

3.968 


.329 

.274 

Hi 

.044 

.30 

3.35 

2.72 

.638 

.084 

3.417 


.337 

.281 

.056 

.008 

.25 

3.18 

2.56 

.640 

.006 

3.132 

.059 

.306 

.252 

.056 

6 * 10“ 4 

.20 

.983 

.328 

.682 

2 * 10" 4 

.855 

.128 

.090 

.031 

.059 

2xl0- 5 

.15 

— 

.920 

.273 

.650 

4 * 10“ 6 

.272 

.649 

.086 

.031 

.056 







DGR = 

1 





0.40 

1.62 

1.59 


2.79 

4.567 





0.127 

.35 

1.63 

1.60 


.674 

2.349 



I 


| ■ 

.30 

1.64 

1.61 

.031 

.084 

1.661 




.003 


.25 

1.62 

1.59 

.031 

.006 

1.588 

.039 



.003 

6xl0" 4 

.20 

.346 

.322 

.027 

2xl0” 4 

.234 

.112 

.024 

.021 

.002 


.15 

.181 

.156 

.024 

4xl0 -6 

-.359 


.019 

.017 

.002 

4xlO' 7 
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TABLE 10.- VAN BAVEL-CR0P-R0G0WSKI STANDARD 10-DAY TEST RESULTS 


Cumulative centimeters for 10 


Eu Drainage 


7.40 

5.18 

2.20 

7.37 

5.15 

2.20 

6.76 

4.51 

2.22 

4.36 

2.01 

2.35 

3.07 

.65 

2.42 

2.86 

.61 

2.25 


Water 

loss 




BALANS 


30 


17.03 

9.28 

6.79 

4.25 

2.81 

1.82 


Centimeters for last day 


E v Drainage 



0.032 

0.727 

0.533 

0.190 

.051 

.720 

.526 

.191 

.091 

.551 

.351 

.199 

.117 

.286 

.073 

.213 

.258 

.282 

.073 

.209 

1.03 

.259 

.070 

.190 



0.40 

5.33 

3.94 

tm 

9.83 

15.16 

0.003 

0.532 

0.410 

.35 

5.34 

3.95 

11111 

1.99 

7.28 

.05 

.531 

.410 

.30 

5.26 

3.88 

1.39 

i 

.12 

5.34 

.05 

.513 

.392 

.25 

3.85 

2.40 

1.44 

.006 

3.75 

.103 

.248 

.123 

.20 

1.97 

.46 

1.52 

2xl0’ 4 

1.79 

.183 

.181 

.050 

.15 

1.85 

.43 

] 

1.43 

7xl0‘ 6 

1.02 

.834 

.170 

.048 



0.283 

.141 

.012 

6xl0” 4 

2xl0" 5 

7xl0" 6 


i 

0.40 

3.35 

2.71 

0.64 

9.98 

13.34 

0.013 

0.337 

0.281 

.35 

3.36 

2.73 

.64 

2.02 

5.364 

.008 

.337 

.281 

.30 

3.35 

2.73 

.64 

.119 

3.46 

.016 

.334 

.281 

.25 

3.17 

2.54 

.64 

.006 

3.12 

.052 

.291 

.236 

.20 

.98 

.30 

.68 

2 * 10‘ 4 

.853 

.131 

.091 

.031 

.15 

1 

.92 

.27 

.65 

1 

O 

i —* 

X 

.266 

.654 

.097 

. 

.031 

i 



0.40 

1.64 

1.61 

.35 

1.64 

1.61 

.30 

1.64 

1.61 

.25 

1.63 

1.60 

.20 

.305 

.281 

.15 

.181 

.156 



1.607 

.101 

-.375 


0.019 

.009 

.010 

.031 

.104 

.556 



0.317 
.158 
.003 .012 

.003 2xl0“ 4 

.002 2xl0 -5 

.002 7xlO" 7 
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Figure 18.- The soil water profile on the 10th day using 
the Rogowski model for soil water characteristics. 
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Solar radiation (mj/nr) 


Figure 19.- The variation of ET for 10 days versus the solar radiation and 
initial soil water profile using Ghosh's model for soil water 
characteristics. 




Soil water (e) 



Figure 20.- The soil water profile on the 10th day using the 
Ghosh model for soil water characteristics. 
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water loss. In addition,, the BALANS values are smaller for the Rogowskl model 
and largest for the regression model. 

3.3 FALLOW SIMULATIONS 

Fallow conditions were also simulated. This was accomplished by deleting from 
the program some of the terms Involving the crop and by setting LAI and RD to 
zero. The resulting daily evaporation for the 10-day period for the different 
Initial conditions using the regression model are shown in figure 21. The 
cumulative evaporation for the 10-day period is shown In figure 22. These 
simulation results Indicate that the three stages noted earlier for ET are also 
present for evaporation alone. A basic difference In the simulated evaporation 
under a good crop cover (LAI « 3.5) and for fallow condition is that under crop 
cover the E v is nearly always constant, while for fallow conditions It starts 
higher and ends lower. 

The changes In the E y tor 10-day values as the solar radiation (DGR) Is 
increased or decreased Is Indicated In figure 23. These results are similar to 
those given for ET shown In figure 6. However, for comparable DGR, the E v for 
10 days is lower. Also, the response differs since the E v approaches a maximum 
value as DGR Increases. The evaporation on day 10 for different DGR and 
0 values Is shown In figure 24. This response differs significantly from that 
for ET which Is illustrated in figure 8. The 10-day soil water profiles are 
presented in figure 25. These profiles look more realistic than do the 
profiles under crops (see figs. 16, 19, and 20). 

3.4 THE BALANS EVALUATION AND RC(I) INCONSISTENCY 

To Investigate the reason the BALANS values increase as the soil becomes dryer, 
the program code was modified to print out the RC( I ) values and CUMRC along with 
the other output information. Inspection of these values for various simula- 
tions indicated that, as the Initial profile was made drier, the absolute value 
of CUMRC became progressively smaller than the value of CUMTR, and it finally 
became positive for the dry regime. This positive value suggests that the crop 
was taking water from the air and putting it in the soil. In addition, as the 
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Figure 21.- The daily E y for the 10-day period for different initial soi 
water profiles for fallow fields and standard atmospheric conditions. 
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lure 22.- The cumulative E y for the 10-day period versus the initial soil 
water profile for fallow fields and standard atmospheric conditions. 






Figure 23.- The changes in the E y for the 10-day period versus the solar 
radiation (DGR) and initial water profile for fallow fields. 
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Figure 24.- The E y on the 10th day versus the initial soil water profile 
and different solar radiation amounts for fallow fields. 
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Figure 25.- The soil water profiles on the 10th day for 
fallow fields using the regression model. 
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soil profile became drier, positive values for RC(I) began appearing near the 
surface and became larger as drier Initial profiles were used. This behavior of 
RC(I) Is Inconsistent with what Is expected to occur, which Is negative RC(1) 
values and absolute CUMRC values always nearly equal to absolute CUMTR values. 

An analysis of the mathematical equations (4), (5), and (6), which determine the 
RC(1) values. Indicates that It Is possible for RC( I > to be positive, especially 
when TRC Is small. The term Involving TRC Is multiplied by SRCR; so. In order 
to test this hypothesis, SRCR was Increased and further simulation performed. 
These runs Indicated that CUMRC and BALANS values were progressively Improved as 
SRCR was made larger. 

In order to further evaluate the nature of the RC(I) Inconsistency, the program 
code was modified so that positive RC(I) values were set to zero, and the model 
was run for the range of 8 values under standard conditions. The simulation 
results showed that, at the wet boundary, CUMRC was nearly equal to CUMTR, and 
the values were nearly the same as those provided by the original code. 

However, as drier Initial conditions were simulated, the BALANS values became 
progressively larger, reaching a maximum value at 9 = .30 of - 21.3 cm with 
CUMRC the greater. As a further check, the relation below was substituted for 
the original code. 

RC(I) = RF( I)*TRC (16) 

As expected, this equation provided values of CUMRC just about equal to CUMTR. 

In addition, the BALANS values became small for all 9. Examples of the model 
response using equation (16) are presented in Appendix B. 

One conclusion that can be drawn from the above results is that the original 
program provides positive values of RC(I) in some layers and negative values in 
other layers. Furthermore, the positive values become progressively greater 
than the negative values for progressively drier conditions with the largest 
positive values near the surface. The Interpretation here Is that the model 
simulates plant root uptake of water in some layers (negative values), and it 
simulates a loss of water to the soil In other layers (positive values). The 
overall result, at a given time. Is that the soil water In the profile near the 
surface is too much when compared to the amount it would be If the RC were In 
balance with the ET (i.e., small BALANS value). 
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4. SUMMARY OF SIMULATION RESULTS 


MATBAL1 Is a computer model that predicts the evapotransplratlon and the soil 
water profile as a function of time. The computer program solves the nonlinear 
partial differential water transport equation numerically using the CSMPIII. 
This latter program Is easy to program and use. 

Although the model Is quite complicated with a number of empirical equations 
and coeffl clents , the output obtained from the sensitivity study appears quite 
reasonable and realistic. A number of the response curves agree with the 
results of empirical studies. 

The sensitivity analysis did Indicate several unrealistic responses In the 
Intermediate and dry regimes In both the ET for 10 days and the BALANS 
values. The cause of these responses were located In the algorithm that 
determines the water uptake by the roots. At low-water amounts, this algorithm 
took water from the soil at deep layers and put water in the soil In the near- 
sv»*face layers; the drier the soil the more pronounced this effect. For 
extremely dry conditions, water was essentially taken from the air and put Into 
the soil. 

Simulations using an algorithm that equated the root uptake to the ET provided 
responses that were similar in the wet regime but more realistic in the Inter- 
mediate and dry regimes. Presented below are the specific responses of the 
model to parameter or atmospheric changes: 

1. The model simulates the diurnal variation in soil moisture, the amplitude 
of which decreases with depth as expected. 

2. In general, the model responses to changes in atmospheric evaporativity 
appear reasonable and realistic. 

3. Both crop and fallow cases reproduce three- stage drying. 

4. Total water loss becomes progressively less than ET as drier initial 
profiles are simulated, an unrealistic occurrence. 
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5. BALANS values become progressively greater as drier Initial conditions are 
simulated, and they eventually become much too large. 

6. The response to Increases In wlndspeed and LAI/RD In the dry regimes do 
not appear realistic. 

7. The response In the wetter regimes to Increases in LAI/RD are small over 
most of the range of values. 

8. The response to percentage changes in SRCR are small. Decreasing the 
value In the drier regime can change the sign of the ET indicating that 
water Is taken from the air and forced into the soil. Increasing the 
value of SRCR decreases the transpiration, but it also decreases the 
BALANS value. 

9. Responses to changes In RLVSWP gives very little change. 

10. Responses to changes of WPCRM are negligible. 

11. Soil water profile changes are negligible In the drier regimes and provide 
unrealistic profiles. 

12. Changes In hydrologic properties were investigated by using three models: 
regression, Ghosh, and Rogowski. 

a. The T, E v , and ET values for the wet and dry soil cases show little 
differences among the models; however, in the intermediate regimes, 
the values vary quite significantly. 

b. Drainage and water loss are least for the regression model and 
greatest for the Rogowski model. The values are large at the wet 
boundary and then become progressively less for increasingly drier 
conditions. Although the drainage becomes insignificant in the 
intermediate and dry regimes in all three models, a large difference 
exists between models in wet regime. 

c. The 10-day soil water profiles are in general similar, but none of 
them show realistic water losses near the surface in the dry regime. 

13. The response of E v for 10 days to changes in solar radiation for fallow 
conditions are similar to the ET changes in the crop case In wet regimes 
but somewhat less In the intermediate and dry regimes. 
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14. The 10-day soil water profiles for the fallow case show large surface 
drying with final surface values of the different profiles very close 
together with steep gradients. 

15. As a function of soil water amount and solar radiation values dally E v 
values for the fallow case have responses similar to those found by 
Denmead and Shaw. 

16. The three-stage drying response Is also provided by the modified model. 

17. The ET for 10-day response In the modified model to changes in DGR as a 
function of 0 are similar In character to the fallow response, but larger 
In value. 

18. The ET for 10-day response to changes In LAI/RD are more realistic than 
the original model. 

19. The 10-day profiles show surface layer drying In the modified model which 
are more realistic than the profiles from the original model. 
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5. CONCLUSIONS AND RECOMMENDATIONS 


The responses of the model have been tested for a range of values for most of 
the atmospheric, crop, and soil parameters. In particular, the response to 
rainfall was not Investigated systematically, but It appears to be realistic In 
‘ general use. Most of the responses to the tests appear to be realistic; 
however. It was determined that the logic that was related to the root uptake 
of soil water did not appear to give reasonable responses In the Intermediate 
and dry regimes. When the logic was modified to relate total root uptake 
directly to transpiration, the model provided more realistic responses. 

In general, the ratio of the percent change In response to percent change In 
Input Is one or less than one. None of the cases Investigated provided an 
unreasonably large percent change in the response. However, LAI, RD, RL, and 
SA cannot be allowed to be zero because they occur In the denominator of a 
mathematical term. 

Because the positive and favorable aspects of the model surpass the negative 
aspects. It Is recommended that It be tested with field data along with the 
other models. In particular, the model language CSMP III Is flexible, easy to 
modify, and easy to use. 
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FUNCTION 


• ••• 

function 


• ••• 

FUNCTION 


WATBAL1 MODEL LISTING 


VOLUMETRIC 
TVSPl * ( 


VOLUMETRIC 
TVSC1 = ( 


VOLUMETRIC 
TVSP? = ( 


0.010 

0.030 

0.050 

0.070 

0.090 

0.110 

0.130 

0.150 

0.170 

y ,a 


6.230 


WATER CONTENT VS. PRESSURE POTENTIAL 
0.000* -0.7000E*V6) * . 

-0.4000t 4 06) * 

-0.7000t*05) • 

-0. Ib50t 4 05) • 

-0.9500E*04) ♦ 

-0.5700t 4 04) * 

•0.346ot *04) * 

-0. I950f04> * 

•0.1 150t*04) . 

-0.64GOE 4 03) 

-0.3fl00£ 4 03) . , 

-0.2250E*03) • . 

-0.1300E*03) . 

0.250* -0 . foil OE ♦ 02 ) * 

0.270* -0.430E 4 02)* 

0.290* -0.2500F *02) • 

C.310* -0.l5u0E*02) . 

0.330* -o. 7root*on 

0.350* -0.3000E»01) 

0.370* -Q.2e00E *01 ) 

0.390* -0. 1200E 4 0 I ) 

0.410. -0.3000E 4 00> 

0 • <*30 « -0.1000E-01) 

1 . 000 * 0 . 00 ) 

WATER CONTtNT VS. HYDRAULIC CONDUCTIVITY IN M/S 


0 . 000 * 
0.020 
0.040 
0 . ObO 
0.080 
0.100 
0 . 120 
0. 140 
0. 16C 
0.160 
0.200 
0.220 
0.240 

8:168 

0.300 

0.320 

0.340 

0.360 

0.380 

0.400 

0.4?0 

1.000 

WATER 

0 . 000 * 

0 . 010 * 

0.030* 

0.050* 

0.070* 

0.090. 

0 . 110 . 

0.130* 

0.150* 

0.170* 

0.190* 

0 . 210 * 


0.0000000 
0. lbbbOOOE- 
0.527800UE- 
0. lfiblOOut- 
0.63B9000E- 
0,227800 OE* 
0.8333000E- 
0.291 7000E- 
0. lOOOOOOt- 
0. 361 1000E- 
0. 1250000t- 
0.**4 *mOOO£” 


) • 

7) 

7) 

6 ) 

6 ) 

S) 

5 ) 

4 ) 

3) 

3) 

2 ) 

2 ) 


0. 1500000E-U ) 

8:tJ»888f:HI 

0.5278000E-10) 

0. 1611000E-09) 

0.5278000E-09) 

0. 1 722000E-08) 

0.5556000E-08) 

0. lV4**OOOt~C7) 

0.11 1 10UOE-07) 

0. 1 11 1000E-07) 

CONTENT VS. PHtSSURE POTENTIAL 
- 0 • aOO OE 4 Ob ) » 

0. 1 OOOE ♦ 06 ) . 

O.SOOOt 4 05) ♦ 

0 . 2«*50t ♦ 04 ) . 

0. 15b0t*04) . 

0 . 9b00t ♦ 0 3 ) » 

0 . 64 OOt ♦ 0 3 > • 

0.41 OOL *03) ♦ 

0 . 2600L 4 0 3 ) * 

0. 1 7 00E 4 0 3 ) * 

0. 1070L*03) * 

0.6bO0E 4 02) ♦ 
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FUNCTION 



=1:1,. 

VOLUMETRIC < wit?8°{0NTENT 0 vi. * 
VSC2 * < 0.008* 0.0000000 

.... - 6000E- 


:|:2tl!fc:i 

-o.l6f 

-0.96( 
-0.561 __ 
•0.3300t*0 


•0. 1900E*Q1 ) 


1200 


60 


00 


o.o?o 

8*^6 

S:n 

0.120 

0.140 

0.160 

0.180 

0.200 

8:158 

0.260 

0.200 

0.300 

8:558 

0.360 

0.300 

0.400 

0.420 

158?° 



WRITE (6* 1 _ - 
1100 FORMATPO 1 

00 40 1*1 *NLL 


0.3194000E- 
G.9306000E- 
O.2770OOUE- 
0.03330006* 
0.2333000E- 
0.70030006- 
0. 

0. 

0.17220006- 

O.4503QOUE- 


0, 

0, 


) » 

• • • 

16) 

t • • • 

5) 

f • • • 

5) 

t • • • 

15) 

t • • • 

4) 

• • • • 

14) 

t • • • 

4) 

t Ml 

13) 

t • • • 

13) 

• • • • 

2) 

• • • • 

12) 

}j 

f • • 6 

t • • • 

0) 

0) 

1 • • • 

t • • • 

t • • • 

09) 

1 • • • 

09) 

t f • • 

00) 

00) 

t • • • 

t • • • 

on 

t • # • 

07) 

t • • • 

07) 



... 


... 


TCOM 


DEPTH 


ITHETA* ) 


40 WRITE(6.i200) 

FORMAT (1H *I2*3F10.5) 


I.TCOM(I) * DEPTH! I) • ITHETA (I) 


HTIM 
STIM 


•TIME/3600. 
MOOtHTF ‘ 


EAMOO(HT IME*24.) 

Y*IMPULS (86400. *06400.) 

IF ( Y.LT .0.5) 60 TO 22 
DNUM * DNUM * 1 
22 CONTINUE 

XONUMaFLOAT (ONUM) 

JONUM * W INPUT < 1 *tiNUM) 

?HE f A < ! * 1 * 1 VOL W ( I ) / T COM ( I ) 

CONO ( I ) *AFGEN ( TVSC1.THETAU) ) 

RI8f 1 1 t :tf8f ? i > lUIStAttf TA ( 1 n 
50 S8 N IiT-i4.NL 


THETA (!) * VOLW ( 1 ) /TCOM < 1 ) 
“3NO ‘ 


C0N0(I)*AF6EN( TVSC2*TH£TA ( I ) ) 

PPOT ( f ) *AFGEN ( TVSP2* THETA ( I ) ) 

HPOT ( I ) *PPOT ( I ) -DEPTH ( I ) 

CONTINUE 
00 60 I s 2*NL 

AVCONO ( 1 > * ( CONO ( I - 1) * T COM ( I - 1) *CONO (I ) • TCOM CD)/... 
( T COM C I — 1 ) *TCOM( I ) f 

CONTINUE 

FLUX(NLL) * CONO (NL) 
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90 


00 90 
FLUX ( ] 


lorn -1 > -HPOT < I ) ) *AVC0N0 (D/OIST (1) 


33 


* I>*<H 

CONTINUE 

BEGIN » W INPUT (9.DNUM) 

END * WlNPUTUOtONUM) 

RFT « W INPUT ( 1 I* QNUM) 

RA|N«0.0 

IF (RFT .EQ.O.O) GO TO 33 
URSLOPa < > / < (ENO-BEGf N) *«2> 

OWSLOPs-yPSLOP 
MQPNT;(BCGIN*6 
HEIGHT. (_ _ _ _ _ 

I F ASTiME.G£.eEGlN.ANO.STlHE.LE,MOPNT)HAlN 
(UPSLOP* (ST IME-BEGIN) >/ JbOOOOO.O 
IF(STIME.GT.MDPNT.AND.STlMe.LE.~ 

(OWSLOP* (STIME-ENO) ) /3600000 .0 


jND)/2.0 

!*0 # R^T)/ (ENO-BEGIN) 




>END) RAIN* 


« • « 


CONTINUE 

ol?*1nput<; 


1 • DNUM > 


OGR/86400.*l.E06«24./DL*PI/2.s436.33«DGR/OL 

GR*436 • 33* w I NPUT ( 3 • DNUM > /DL*S I N ( (STIME-12.*DL/2. ) • . • 


44 


500 


501 


•3.141/DL) 
(GR.LE.O.O) 
(HTJME.LE.12, 


GO 
) 


TO 160 
SA=WINPUT(B.DNUM) 


!f(HTIME. 

IF(HTIME.LE.12.) GO T C 44 
IF(STlME.LE.12.)SA=Wu\PUT (R,0NUM-l>*(STIM£*12.)/24.* ... 
(BlNPUT(d.DNUM)-*lNPUT(8.DNuM-l) ) 

IF(STIME,LE.12.) GO TO 44 

SA=tf INPUT (8*DN0M; ♦ (ST IME- 12. > /24.* ( RlNPUT (8*DNUM* 1 ) . . 
WINPUT(8.0NUM) ) 

CONTINUE 

RA = (AL06(2.0/Z0)**2.0)/(0.16*SA> 

OPMAX=tt INPUT (6. DNUM) 

DPMIN=<rf INPUT <7. DNUM) 

DPTC=DPMIN* (DPMAX-OPMIN)*(STlME-5.)/10. 

IF(5TIME.GT.15. )0PMIN=W1NPUT (7*DNUM*1) 

IF(STImE.GT.15.)0PTC=0PMaX-(DPMAX-0PMIN)*(STIME-15.>/14. 

IF(STlME.LT.5.AND.DNUM.oE.2.)DPMAXstfINPUT(6*DNUM-l) 

IF (ST IME. 1 T.5.)DPTC=DPMAA-(DPMAX-DPMIN)*(STIME*9.»/l4. 

HA = 1 . 323 *EXP ( 1 7 . 27 * 0 P TC/ ( 237 . 3 * 0 PTC) > / ( 273 . 16 ♦OPTC) 
TAMAXs. INPUT (4. DNUM) 

TAMIN=y INPUT (b.ONUM) 

TAC=TAMIN* (TAMAX-TAM1N)*(STIME-S.)/10. 

IF (ST IME. GT. 15. ) TAMIN=WINPUT (5.DNUM*1> 

IF (STIME.GT . 15* ) T AC = TAM A a- ( TAM AX* TAM IN) * (STIME-15. ) /l4# 
IF ( ST IME. LT. 5. ANO. DNUM. Gt. 2. ) TAMAX=*1NPUT (4.DNUM-1 ) 
IF(STIME.LT.5.)TAC=TAMAA-(TAMAA-TAMIN>*(STIME»9.)/14. 

T AK = T AC ♦273. lb 

SH= (1154.8*303. 16) /(TAK) 

SKL=(SIGMA*TAA**4) *(0.605*0. 039*SORT( 1410. *HA) ) 
RD=WlNPUT(13fDNUM) 

00 500 1=1. NL 

RF( P =2.0® ( 1/Pi)-0EPTH( I ) /R0**2)*TC0M(I) 

IF (PF < I ) .Ll.O.O) RF< I)s0.0 
CONTINUE 
WPSEFF=0. 0 
DO 501 1=1. NL 

WPSEFF=*(PSEFF*PPOT ( I ) *RF ( I ) 

CONTINUE 

LAl=WlNPUT <12. DNUM) 

ALBC=0. 1240-0. 009988*LA 1*0. 007142*LAI**2-0.000583*LAl**3 

A8SC=0. 0032* 0.3084*LAI-0.05323*LAI **2*0. 003667*LAI**3 
FTSR=0. 9842-0. 6755*LAI *0. l595*LAl**2-0. 0 124 1*LA1**3 
WPOTCR=-WPOTCR 
RL*AF6EN (RLVSttP* WPOTCR) 
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QUALITY 





(273* 16+ TO 

LM«2.49463E06-2.247E03*TL 
LTR s(HL-HA)*LH/RC* 

SHCA*-NRBC*LTR 

ftl*tac-$hca*ra/sh 

TRC*(HL-HA)/<RC«* 1000.0) 

NRBS*GR* ( 1 . O-ALBC-AdSC ) *F TSR* ( SKL-LBRC > 
FPS*0.92l-0.002o2*rL*0.p0308*TL**2 
CVS* (EPS/ (EPS* 1.0)1 *NRB5/ (LH* 1000*0 
EVS*EVS*EXP (PPOT ( 1 ) / (46.97* T AK) ) 
IF(TRC.LT.O.O) TRCaO.O 
IF (EVS.LT.O.O) EVSsO.O 
WPOTCRsWPSEFF»dPC«MN-TRC*SRCR/LAl 
00 502 I*1.NL 

RC(I)*(WPOTCR-WPCRMN-PPOT (I) > *RF ( 1 ) *LAI/SRCR 
502 CONTINUE 

EVTRsEVSOHC 
GO TO 170 

160 CONTINUE 
EVS=0.0 
TRCsO.O 
EVTR*0.0 

00 161 1*1. NL 
RC< I)*0.0 

161 CONTINUE 
170 CONTINUE 

(DETAIN r INTGRL (U.O, RAIN-INFILT) 

INCAP = (O.-HPOT(l) )*0.5*(SATCON«COND(1> > 

IF (RAIN. GT. 0.0) GO TO 5b 
IF (DETAIN. LE. 0.0) GO TO 66 
INF ILT* INCAr 
GO TO 77 
CONTINUE 
DETAIN = 0.0 
INFILT-O.O 
GO TO 77 
CONTINUE 
INFILT = 

IF (RAIN, 

CONTINUt 

FIUX(1)=INFILT-EVS 
00 100 I = l.NL 

NFLUX ( I ) sFLUX ( 1 ) -FLUX ( I ♦ 1 ) *RC ( I ) 

CONTINUE 

V0L4= INTGRL ( I VOL W. NFLUX. 18) 

COMRN s INTGRL (0.0 . RAIN) 

CUMINF * INTGRL ( 0.0 . INFILT ) 

CUMEV = INTGRL ( 0.0 . EVS ) 

CUMTR = I,NTgRL(O.O.TRC) 

CUHETRs INT6RL (O.O.EVTH) 

CUMORN = INTGRL ( 0.0 . F LUX (NLL ) ) 

CUH*TR = 0.0 
00 110 1=1. NL 
CUMrfT«= CUHrtTR ♦ VOLW(I) 

110 CONTINUE 

ZBHJS=IMPULS( 66400. . 86400.) 

IF (ZRHJS.LT. 0.5) GO TO 86 
INF(DNUM-1)=CU*INF 

RN(DNUM-1)*CUMRN 
EV (0NUM-1)=CUHEV 
TP (DNUM-i ) sCUMT R 
ETR(DNuM-l) *CUMtTR 


/ OIST(l) 


66 


55 


77 


100 


INCAP 

LT. INCAP. 


ANO. DETAIN. LE. 0.) INFILT*RAIN 
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s 

% 


0«AIN(DNUM-n «CUMUPN 
DINF«CUMlNF-INf (0NJM-2) 
OwNsCUHRN-kN ( ON JM-2 1 
OEV*CUM£V-EV (UNJM-2 ) 
OT«*CUMTR-TR t ONOM-2 ) 
0£TR*COMET«-ETR OJNJM-2) 
OORNsCUMORN-OR A IN ( ONUM-a J 


88 CONTINUE . „ , 

8ALANS 3 CUHdTw - HATER - 
Z*ImPulS( 10800.0. I OSOO . 0 » 
IF(Z.LT.O.S) 60 TO 99 
222 CONTINUE 


CUMINF ♦ CUMETR ♦ CUMDRN 


1300 P 8 A 5 S n * 1 38 fi l A « 1 EfSV T * I i PI : 0 f 9 N V?Af T i M 5 : ?*r 0 . 1 . 

• XDNUM s «.FU.0.» STIME s «*F7.4*« T L 3 «.F4,1> 
KRITE(8.1400) 

1400 FORhaTi»0 I ».5A**iJEPTM».10x.*THETA».11a.»PPOT*. 11A. 
•FLUX* *9X* » NET FLUX* • 10X. *POuTF • . 10X. •POOTUPTA*E • ) 

00 150 1*1 *NLL 

150 tffll TE (6. 1500) l.OEPTH(I) .THFTA(I) .PPCT 1 1) .FuUX ( I) . 
NFLUX(I) *HF ( 1 1 ,PC<I) 

1500 FORMAT { I3.7E15.4) 

*RITE <6* 1600) 

1600 FOP 4AT ( 1H ////) 

99 CONTINUE 


TIMER FINTlMa 8640 00 . o P«0ELs 86400 .. OUTOEL=720 0 . 0 *GtLTslO. 
PPINT JGNuM* XDNoM. QPN. OInF. uE V .OTP *D£TP . OQRNt CUH<iT9*« # . 
8AL ANS * CUMRN. CUMOPN. CUM I nF .CUMEV .... 

CUMTW.CumETN.DELT 


PAGE SHAOE * (0.15*0.35) 

OUTPUT TmE T A ( 1 3 1 .THFTAQ2) • Th£T A ( 1 1 ) » TriE T A ( 1 0 ) *Th£T A (9) *••• 

theta (8) .theta ( 6i .theta t ai 
MPTHOO stiff 

E 40 


♦ <EaTH£W 

Input 

data 

• STOkEj 

In ARRAY *InPuT < 13*37) 


• JNM 

input 

OL 

OGR 

Th a * 

TMIN 

OMAX 

OHIN 

SA 

BEGIN 

ENG 

RFT 

LAI 

121 . 

14 .4 

20.0 

25.0 

10.0 

10.0 

6.0 

3.0 

0.0 

0.0 

0.0 

3.0 

122. 

I 4 • 4 

20.0 

25.0 

10.0 

10.0 

e.o 

3.3 

0.0 

0.0 

0.0 

3.1 

12**. 

1 4 .5 

20.0 

25.0 

10.0 

10.0 

8.0 

3.0 

0.0 

00.0 

00.0 

3.2 

124. 

N.S 

20.0 

25. 0 

1C.0 

10.0 

e.o 

3.0 

0.0 

0.3 

0.0 

3.3 

125. 

14.5 

20.0 

25.0 

10.0 

13.0 

3.0 

3.3 

C • 3 

0.3 

0.3 

3.4 

126. 

14.5 

20.0 

25.0 

10.0 

1 0 . 0 

c. 0 

3. 3 

0.3 

0.3 

0.3 

3.4 

127. 

W.5 

20.0 

25.0 

10.0 

10.0 

*. 0 

3. 3 

0.3 

J . 3 

J . 3 

3.3 

126. 

i 4 , 6 

20.0 

25.0 

10.0 

10.0 

e.o 

3.3 

0.3 

0.0 

0.0 

3.6 

120. 

U.6 

20.0 

25.0 

10.0 

1 0 . 0 

8.0 

3.0 

0.3 

0.0 

0.0 

3. 7 

130. 

0. 

14.6 

14.6 
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25.0 

10.0 

10.0 

10.0 

10.0 

6.0 
8 . 0 
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3.3 

0.3 

0.0 

0.0 
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APPENIDX B 
MODIFIED WATBAll 


This appendix presents the results of further simulations that were made In 
order to better evaluate the effects of the root uptake function Inconsistency 
on the previous simulation results. These simulations were made using the 
equation: 

RC(I) • RF(I)*TRC (B-1) 

which equates the root uptake with the transpiration. The simulation results 
are presented In figures B-l through B-7. The daily values of T, E v , and TR 
are shown In figures B-l through B-3. These figures correspond to figures 2 
through 4 for the original model. Comparing these figures indicates that the 
results are similar and that the two drying stages are evident in the new 
simulations. However, in the latter curves, the falling stage commences sooner 
and the drop In ET and T are quicker and larger. Another difference Is that E v 
eventually falls below the T value in the modified rersults. The cumulative ET 
response curves are presented in figure B-4. 

The change of ET for 10 days with the daily radiation amount is shown in figure 
B-5 which corresponds to fiqure 7. Comparing the results of the two models 
shows that the modified model has a different response- in the intermediate and 
dry moisture regimes. In particular, the dry boundary does not increase as 
rapidly with increasing DGR values. In the intermediate zone, the ET for days 
increases up to a certain value and then remains more or less constant above 
that DGR value. The humps on the curves for 0 = .30, .35, and .40 appear to 
reflect the rapid falloff from the constant stage. Increases of DGR beyond 30 
mJ/m^ should result in the ET for 10 days eventually increasing aqain. This is 
suggested by the curves in figures B-6 and B-7. 

The manner in which ET is divided into T and E v is shown in figures B-5 and B-6 
(note change of scale in figure B-6). The response curves In these figures 
Indicate that the E v Increases regularly with increasing DGR. T increases up 
to a certain value of DGR, and then it increases In DGR, apparently reflecting 
the change from the constant to the falling transpiration stage. The Indica- 
tions are that it should Increase again for further Increases of DGR (see 
figure B-7). The evapotransplratlon on the 10th day, as related to DGR 
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conditions usiru) t ht' inod i f iod V.in Ravol model 
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il water conditions using the modified Van Bavel model. 
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F in U re B-4 . - The cumulative ET versus the initial soil water 
profile usinq the modified Van Bavel model. 
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Figure B-7.- The variation of the and T f o r 10 days versus solar 

radiation for = .25 and .15 using the modified Van Bavel "iodel . 


and 0, Is Illustrated In figure B-8. lhe curves are similar to those presented 
by Denmead and Shaw (ref. 12) which were obtained from experiments with corn 
grown in larqe pots. 

The ET for 10 days for different LAI/RD values are shown in figure B-9. 
Comparing these results with those in flqure 12 indicates that the dry regime 
response is now more logical . The intermediate reqlme has also changed, and 
proqressi vely higher ET for 10-day values occur as the LAI/RD values 
increase. How these values are divided into T and E v are shown in figure 
B-10. As would be expected, T increases with increasing LA1/R0, while E v 
decreases. 

The final comparison is in the 10-day profiles. These are shown, for the 
regression model, in figures 16 and B-ll. The model profiles in fiaure 8-11 
are more realistic than those presented in figure 16. These latter profiles 
reflect the fact that water is apparently simulated as being taken from lower 
soil layers and transferred to surface layers by the root system. Comparing 
fiqure B - 1 2 for the Ghosh model with the similar conditions in fiqure 20 
indicates that the same types of changes are indicated. 

These comparisons of the simulation results of the modified Van Ravel model 
with the results of the oriqinal model show that some significant differences 
are indicated in the output. These differences become greater as drier reqimes 
are simulated, and they generally appear to reflect the indications that water 
in the original model is simulated as beinq extracted from the soil by the 
lower parts of the root system and returned to the soil in the near surface 
layers. 

In addition, the results from the modified Van Ravel model appear more 
realistic in the drier reqimes than the original model when the experimental 
field data are considered. 
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Soil water (9) 



Figure B- 12.- The soil water profile characteristics on the 10th day 
using the modified Van Bavel model and Ghosh's model for soil water. 
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